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The 't Hooft six quark flavor mixing interaction (TV/ = 3) is bosonized by the path integral 
method. The considered complete Lagrangian is constructed on the basis of the combined 't Hooft 
and f/(3) x (7(3) extended chiral four fermion Nambu - Jona-Lasinio interactions. The method of 
the steepest descents is used to derive the effective mesonic Lagrangian. Additionally to the known 
lowest order stationary phase (SP) result of Reinhardt and Alkofer we obtain the contribution from 
the small quantum fluctuations of bosonic configurations around their stationary phase trajectories. 
It affects the vacuum state of hadrons at low energies: whereas without the inclusion of quantum 
' fluctuations the vacuum is uniquely defined for a fixed set of the model parameters, fluctuations give 

■ rise to multivalued solutions of the gap equations, marked at instances by drastic changes in the 

' quark condensates. We derive the new gap equations and analyse them in comparison with known 

. results. We classify the solutions according to the number of extrema they may accomodate. We 

' find up to four solutions in the < m„,ms < 3 GeV region. 
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I. INTRODUCTION 
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00 ' The global ?7l(3) x Ub.{^) chiral symmetry of the QCD Lagrangian (for massless quarks) is broken by the Ua{^) 
Adler-Bell-Jackiw anomaly of the 5*^7(3) singlet axial current 97/^75(7. Through the study of instantons it has 

been realized that this anomaly has physical effects with the result that the theory contains neither a conserved U{\) 
' quantum number, nor an extra Goldstone boson. Instead the effective ^Nf quark interactions arise, which are known 
. as 't Hooft interactions. In the case of two flavors they are four-fermion interactions, and the resulting low-energy 
theory resembles the old Nambu - Jona-Lasinio model . In the case of three flavors they are six- fermion interactions 
which are responsible for the correct description of rj and 77' physics, and additionally lead to the OZI-violating effects 



I 
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^ • -Cdet = ^(detgPflg + detqPLq) (1) 

^ ' where the matrices Pr,l = (1 i 75)/2 are projectors and determinant is over flavor indices, 
•'-j . The physical degrees of freedom of QCD at low-energies are mesons. The bosonization of the effective quark 
/\ [ interaction (|l|) by the path integral approach has been considered in where the lowest order stationary phase 
^ approximation (SPA) has been used to estimate the leading contribution from the 't Hooft determinant. In this 
■ - - ' approximation the functional integral is dominated by the stationary trajectories rg^(a;), determined by the extremum 
condition 5S{r) = of the action S{r). The lowest order SPA corresponds to the case in which the integrals 
associated with 6'^S{r), for the path rgf^{x) are neglected and only S{rgf^) contributes to the generating functional. 
The next natural step in this scenario is to complete the semiclassical result of Reinhardt and Alkofer by including 
the contribution from the integrals associated with the second functional derivative S'^S{ri,f^), and this is the subject 
of our paper. 

An alternative method of bosonizing the 't Hooft determinant has been reviewed in . The special path integral 
representation for the quark determinant has been obtained by considering Nc as an algebraically large parameter. 
One should not forget that the 't Hooft 's determinant interactions are induced by instantons and only can be written 
in the simple determinantal form (^ in the limit of large number of colours - otherwise the many-fermion interactions 
have a more complicated structure. For our calculations it means, in particular, that the terms of the bosonized 
Lagrangian induced by the 't Hooft 's determinant interaction (|l|) should be of order 1, corresponding to the standard 
rules of Nc counting. In this respect it is worthwhile to note that we have found that the meson vertices, induced by 
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the 't Hooft's determinant interactions in the lowest order SPA and the term from the integral associated with the 
second functional derivative S'^S{rgf^) have the same iVc-order and should be considered on the same footing. 



II. PATH INTEGRAL BOSONIZATION 



To be definite, let us consider the theory of the quark fields in four dimensional Minkowski space, with dynamics 
defined by the Lagrangian density 

^ = %JL + ^det- (2) 

The first term here is the extended version of the Nambu - Jona-Lasinio (NJL) Lagrangian = Co + C^j^^, 

consisting of the free field part 

Co = q{i-f''d^ - rh)q, (3) 
and the U{3)l x U{3)r chiral symmetric four-quark interaction 

Ant = f [(9Aa<z)' + {qil5Kqf] ■ (4) 

We assume that quark fields have color and flavor indices running through the set i = 1,2,3; Xa are the standard 
C/(3) Gell-Mann matrices with a = 0, 1, . . . , 8. The current quark mass, m, is a nondegenerate diagonal matrix with 
elements diag(7fi„, m^, m^), it explicitly breaks the global chiral U{3)l x U{3)r symmetry of the ^Cji^j^ Lagrangian. 
The second term in (||) is given by (|l|). Letting 

Sa = -qXaq, Pa = qijs^aq, S = SaXa, P = PaXa (5) 

yields 

-^det = + + ~ (6) 

with determinants written in terms of the mesonic type quark bilinears. This identity is a first step to the bosonization 
of the theory with Lagrangian 

The dynamics of the system is described by the vacuum transition amplitude in the form of the path integral 

Z = J VqVqcxp (^i J d^xCj . (7) 

By means of a simple trick, suggested by Reinhardt and Alkofer, it is easy to write down this amplitude as 



Z = J VqVqDaaV(j)aVrlVr'^ exp J d^xC 



(8) 



with 



£' = q{iY'd^ - m - fj + i-f5(t))q - ^ [(ffa)^ + {(paf] + r?(aa + GqXaq) 

+ r^{(j)a - GqijsXaq) - [det(cr + ic/)) + det(cr - i(f>)] . (9) 

Eq.(||) defines the same expression as Eq.(^. To see this one has to integrate first over auxiliary fields rf, rj. It 
leads to 5-functionals which can be integrated out by taking integrals over ctq , and (j)a , and which bring us back to the 
expression (^. From the other side, it is easy to rewrite Eq.(^ in a form appropriate to finish the bosonization, i.e., 
to calculate the integrals over quark fields and integrate out from Z the unphysical part of the auxiliary , scalar 
fields. Indeed, introducing new variables a a + Gri, cf) + Gr2, and after that ri 2ri — cf/G, r2 2r2 — (f'/G 
we have 



Z^ J 'D(Ta'D(t)a'Dq'Dqexp ^ J d'^xCq{q,q,a,(t))^ J 'Dria'Dr2aexp Jd'^xCr{(J,t 



4>,n,r2)] (10) 



2 



where 



Cq^ q{i"f^d^^- m- a + i^5(j))q, (11) 
Cr = 2G [{riaf + {r2af] - 2{riaOa + Tsa^a) " ^ [det(ri + ir2) + clet(ri - ir2)] . (12) 

The Fermi fields enter the action biUnearly, we can always integrate over them, because in this case we deal with the 
standard Gaussian type integral. At this stage one should also shift the scalar fields (Ja ^ cFa + by demanding 
that the vacuum expectation values of the shifted fields vanish < 0|(Tq|0 >= 0. In other words, all tadpole graphs in 
the end should sum to zero, giving us gap equations to fix parameters A^. Here — to^ — tti^, with denoting 
the constituent quark masses |]. To evaluate path integrals over ri 2 one has to use the method of stationary phase, 
or, after the formal analytic continuation in the time coordinate — ixo, the method of steepest descents. Let us 
consider this task in some detail. 

The Euclidean (imaginary time) version of the path integral under consideration is 

J{a,(f>)^ J 'Dria'Dr2aeyip (^J d'^x£ricr,4>,ri,r2)^ ■ (13) 

This integral is hopelessly divergent even if k = 0. One should say at this point that we are not really interested in 
(|l|) but only in its analytic continuation. Let us suppose we analytically change £^(0", , r2) in some way such that 
we go from this situation back to the one of interest. To keep the integral convergent, we must distort the contour 
of integration into the complex plane following the standard procedure of the method of the steepest descents. This 
method gives the first term in an asymptotic expansion of J(cr, </>), valid for h —^ 0. We lead the contour along the 
straight line which is parallel to the imaginary axis and crosses the real axis at the saddle point r^^ . It is in the sense 
of this continuation that the integral J(cr, 0) of ( p^ ) is to be interpreted as 

J{a,(j))= / 'Dria'Dr2a exp [ d'^xCr{<7,(t>,ri,r2)] ■ (14) 



St 



Near the saddle point r"^. 



Cr^Cr{rst) + l^rcC'^^{rg^)f(3 (15) 



2 

a, (3 



where the saddle point, r"^, is a solution of the equations £^(ri,r2) — determining a flat spot of the surface 
Cr{ri,r2). 



2Grf - (a + A), - f ^a6c(r?rj - r^r^) = 

2Gr^-ct)a + ^AabaT\r'2 = 



(16) 



This system is well-known from The totally symmetric constants, Aabc, come from the definition of the flavor 
determinant: detr = Aabc'r°"r^ r'^ , and equal to 



Aabc — -^^ijk^mnl{^a)im{\)jn{^c)kl- (17) 



They are closely related with the i7(3) constants dabc- We use in ( |15D symbols f° for the differences (r° — r'^^). To 
deal with the multitude of integrals in ( [1^ ) we define a column vector f with eighteen components = (^17^2) ^^d 



with the matrix C'^Ar^^) being equal to 



O(^st) = 4GQ.„ Q., = ( 3« f f . ) ■ (1^) 



Eq.(O) can now be concisely written as 



^The shift by the current quark mass is needed to hit the correct vacuum state, see e.g. 
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J(a, 0) = exp (^y" d-'xCrirst)^ ^ exp (^G J d''xf^Q{rsi)f^ [1 + 0{n)] . (19) 

Our next task is to evaluate the integrals over fa- Before we do this, though, some comments should be made about 
what we have done so far: 

(1) The first exponential factor in Eq.(p^) is not new. It has been obtained by Reinhardt and Alkofer in A bit 
of manipulation with expressions ([l2|) and (jlq) leads us to the result 



C(rst) = ^ {G[{rl,,r + (r^3t)2] - 2[{a + A)„r^,t + d^^r^^^,]} . (20) 



One can try to solve Eqs.(p^) looking for solutions r°g|_ and r^g^ in the form of increasing powers in (Ta, 4'a 

^tst ^ha + h^^H a, + + h^;^lM, + . . . (21) 

'^2st = /^l'V6 + /iitVc + ... (22) 



K, h^2, and h^^^ 



Putting these expansions in Eqs.(h6f) one can obtain the series of selfconsistent equations to determine the constants 



SG 

The other constants can be obtained from these ones, for instance, we have 



2GK - Aa - ^AabchbK = 0, (23) 

2G (^6ac - ^Aacbhl^ hg^ = 6ae, (24) 
2G (Sac + ^Aacbhb) hf} = 5ae. (25) 



, (1) _ (1), (1). (1) . _ . (2) _ _^l,Wf,(2)ui2)A _ ;,(3) _ _^/,(2);,(2)> (1) . _ 

'afcc ~ g "'aa"'iyi, "-cc -^abci ""abc ~ g "-aa "-fji, "-cc ^abcJ "-abc ~ ^ "•ao "-{jj "•cc ^abc- 



As a result the effective Lagrangian (^0|) can be expanded in powers of meson fields. Such an expansion (up to the 
terms which are cubic in a a, 4'a) looks like 

Crirst) = ~2haaa - h[\^cTaab - hf^Ua^^b + 0(field3). (27) 

(2) Our result ( |l9|) has been based on the assumption that all eigenvalues of matrix Q are positive. It is true, for 
instance, if k = 0. It may happen, however, that some eigenvalues of Q are negative for some range of parameters 
G and k. In these cases there are no conceptual difficulties, for from the very beginning we deal with well defined 
Gaussian integrals and the integration over the corresponding fg simply does not require analytic continuation. 

We now turn to the evaluation of the path integral in Eq.(|l9|). In order to define the measure T^fa more accurately 
let us expand in a Fourier series 

oo 

fa{x) = ^Cn,a^n{x), (28) 
n=l 

assuming that suitable boundary conditions are imposed. The set of the real functions {ipn{x)} form an orthonormal 
and complete sequence 

oo 

d'^xipn{x)ipm{x) ^ 6„m, ^ (p„{x)(p„{y) ^ S{x - y). (29) 



n=l 



Therefore 



J I?f„exp (^2G J d'^xr^Q{rst)r^ = J rfc„,„ exp {2G ^ c„,cA^^c™,/3} 



G 



/det(2GA 

The normalization constant G is not important for the following. The matrix A"^ is equal to 



a/3 \ 
nm I 



(30) 
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KL^ J d^X(pn{x)Qaf:i{x)tpm{x). (31) 



From (|18D and (|29|) it follows that 

■ ,(i)-i 



2GXf^=l^'"'^ (^5,p5^^^ + j d^x^r.{x)F,p{x)^m{x)^ (32) 



with 



Only the matrix F„p depends here on fields a, (p. By absorbing in C the irrelevant field independent part of 2GA"^, 
and expanding the logarithm in the representation det(l + F) = cxp tr ln(l + F), one can obtain finally for the integral 
in dll) 

J(a, 0) = C'e^'-, Sr ^ d^xl Lr{r^^) + - ^ ^^tr [F„"^(rst)] ^m{x)^M ■ (34) 
J I — 1 I 



The sum over m in this expression, however, is not well defined and needs to be regularized. One can regularize it by 
introducing a Gaussian cutoff M damping the contributions from the large momenta k'^ 

X: vMVmix) = <5(0) ^ ^ exp [-— j = ^ (35) 

This procedure does not decrease the predictability of the model, for anyway one has to regularize the quark loop 
contributions in (|l^). Alternatively, following ideas presented in [||, one can introduce the Ansatz: 

Sr^j d\ \cr{r,t) + ^ E [F^pirsi)] | (36) 

proposing that the undetermined dimensionless constant a will be fixed by confronting the model with experiment 
afterwards. 



III. THE GROUND STATE 



Let us study the ground state of the model under consideration, then properties of the excitations will follow 
naturally. To make further progress let us note that Eqs.(p3|) have non-trivial solutions for /iq, /13, /isj corresponding to 
the spontaneous breaking of chiral symmetry in the physical vacuum state with order parameters 7^ {i — d, s). 
We may then use this fact to rewrite Eqs.(|2^) as a system of only three equations to give hi 

2Ghi - = '^t.jkhjhk (37) 

o 

where the totally symmetric coefficients Ujk are equal to zero except for tuds — 1. They are related to coefficients 
^abc by the embedding formula SuJiaAabc^bjeck = Ujk where matrices ojia, and Cai are defined as follows 
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FIG. 1. The couplings G (left panel) and k (right panel) as functions of m„ at fixed values of rria = 572 MeV and other 
parameters in Eq.([4l|). These curves show the typical m„-dependence of the functions if one neglects the fluctuations term in 
Eq.(||). 

^ ( ^/2 V2 ^/2\ ^ / x/2 VS 1 

= \/3 -V3 , oj,a = ^\ V2 -VS 1 | . (38) 

1 1 -2 / V2 -2 

Here the index a runs a = 0, 3, 8 (for the other values of a the corresponding matrix elements are equal to zero). We 
have also ha — Raihi, and hi = coiaha- Similar relations can be obtained for and Aa- In accordance with these 
notations we will use, for instance, that — ujiahi]} . 

A tadpole graphs calculation gives for the gap equations the following result 

^h. + g {h^^ - /.l;^) Aa,.h^ = ^,m.Mml) (39) 



where the left hand side is the contribution from (jS^) and the right hand side is the contribution of the quark loop 
from (|l^) with a regularized quadratically divergent integral Jo(to^) being defined as 

/"OC 7 1 

Mm')=J ^e-'^"p{t,A'), p(t,A2) = l-(l + tA2)exp(-tA2). (40) 




FIG. 2. The couplings G (left panel) and k (right panel) as functions of for the same (as in Fig.l) values of other fixed 
parameters and for the case when the fiuctuations term in Eq.(|3^) is taken into account. 
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FIG. 3. Different classes of m„ solutions of the gap equations, including fluctuations, at fixed G, k values. Four branches of 
solutions. A, B, C and a — b stretch have negative k. Three branches, D, E and along c — d arm correspond to positive k. Bold 
dashes indicate that only one solution exists in this region of G, k. See further details in text. 

The second term on the left hand side of Eq.(|39|) is the correction resulting from the Gaussian integrals of the 
steepest descent method, comprising the effects of small fluctuations around the stationary path, ff one puts for a 
moment a = in Eq.(^), and combines the result with Eqs.(^7|), one finds gap equations which are very similar to the 
ones obtained in [Q (see equation (2.12) therein). At fixed input parameters G, k,A of the model, the gap equations 
can be solved giving us the constituent quark masses rrii as functions of the current quark masses, = mi{mj). 
Alternatively, by fixing rhi, one can obtain from the gap equations the non-trivial solutions rrii = mi{G, A). In 
particular, when rhu — rhd, these equations can be solved for G and k, giving expressions 

/ 27r2 \ m„A^ Jo(m^) - m^A^ Jo(m^) _ / Stt^ \ ^ to^A„ Jo(m^) - TO^AsJo(mg) 

In Figj^ we plot the curves of G and n versus m„ keeping constant A = 0.87 GeV, TOs — 572 MeV, iris — 
200 MeV, m„ = 6 MeV. One can readily see that at given values for {mu,mg) the curves yield unique values 
of (G, k), i.e. the vacuum state is well-defined in this case. 




FIG. 4. The same as in figure 3 for the solutions. The dashes are the constant rus corresponding to the abed curve of 
Figs. 2 and 3. 

Let us consider now the general case which we have when a > in Eq.(|3^). To illustrate the qualitative difference 
with the previous case we put for definiteness a = 1 and look again for the solutions G — G{mu) and k = K(m„) 
with the same set of fixed parameters. The corresponding curves are plotted in Fig.^. If the mass m„ is sufficiently 
low that m„ < mi'"™\ (in the figures denoted by the region left the turning point c) or sufficiently high that 
^(max) ^ ^ ^j^^ figures denoted by the region right to the turning point 6), then there exists again a 



7 



single solution with unique values of (G, k). However, there is now a region for m„ in which m„ " < m„ < niu"'^ , 
where three values of couplings (G, n) are possible. 

Conversely, one can study the solutions: m„ — •m^i^{G^K), rrig = ms{G,K at fixed values of input parameters: 
A, m„ — rhd, rhs). As starting input values for the couplings G and k in the gap equations we take the ones already 
determined along the path abed shown in Fig.|^, obtained at a constant value of the strange quark mass, iris — 572 
MeV. For a chosen value of the set (G, k) we then search for further solutions (m„,TOs) of the Eq.(^9|), displaying 
results in FigJ^ for niu and in Fig.^ for nis, correspondingly The dashed curves are the repetition of the solutions 
encountered in Figj^. The bold dashes in Fig.|^ indicate that we only find this one solution at fixed values of G, n. 
Combining the information of Figs.|| and |[ one sees that one has up to four solutions at fixed G, k. Indeed, travelling 
along the original path abed one observes the following: the branch ab is accompanied by three other branches, marked 
as A, B and C which belong to the same class of solutions. One sees that these solutions have negative k values. From 
the turning point b until the maximum value of G (and corresponding k) (bold dashes) we have no other solutions to 
the gap equations, as already stated. From this maximum G value to the turning point c and further along the cd 
arm we encounter further two branches, denoted by D and E to the solutions of Eqs.(|3^) with same G, k. They are 
positive K solutions. 

This very rich structure of the vacuum solutions implies the possibility of having several different values of the quark 
condensates for the same G, k parameters, embracing also the possibility which has been considered in connection 
with generalized chiral perturbation theory (see e.g. section 4 in We give here only a few examples. At 

G = 4.54 GeV"^,K = 153.04 GeV"^ we have three solutions. The solution m„ = 346 MeV, = 572 MeV, on 
the ed arm has the quark condensates < uu >^/^= —236.8 MeV, < ss >i/-^= —183.5 MeV and the ratio R = 

1 /3 

{< ss > / < uu >) ' — 0.775; the second solution, on the E branch with iriu — 535 MeV, = 732 MeV has 
condensates < uu >^/'^^ -249 MeV, < ss >^/^= -183 MeV and the ratio R = 0.738; the third solution, located at 
the D branch with m„ = 46 MeV^nis = 418 MeV has condensates < uu >i/3= _i3i MeV, < ss >^/'^= -172 MeV 
and R = 1.313. We chose the next example at G = 8.126 GeV~^, n — —544.81 GeV^^, where there are four solutions. 
The solution m„ = 316 MeV, = 572 MeV, on the ab arm has the quark condensates < uu >^/^= —233 MeV, 

< ss >i/3= -184 MeV and R = 0.787; the second solution, on the B branch with m„ = 624 MeV,™^ = 205 MeV 
has condensates < uu >i/3= -249 MeV, < ss >'^/^= -56.7 MeV and R = 0.227; the third solution, located at the 
A branch with to„ = 421 MeV, = 1.84 GeV has condensates < uu >i/3= -14.4 MeV, < ss >i/3^ _85 MeV 
and R = 0.353. The fourth solution, on the G branch, with rriu — 1.394 GeV, to^ = 1.594 GeV has condensates 

< uu >^/3= -230 MeV, < ss >^/^= -120 MeV and the ratio R = 0.524. As a final example we take the solutions 
at G = 11.96 GeV"^ ~ G„iax,K = 371.491 GeV"^ where the branches D and E emerge and are very close to 
each other with niu — 1.56 GeV, — 1.72 GeV. The corresponding condensates are < uu >^/'^= —224 MeV, 

< ss >i/3= -105 MeV and R = 0.467. The other solution is at the path be with to„ = 290 MeV, = 572 MeV 
with condensates < uu >^/^= -229 MeV, < ss >^/^= -184 MeV and R = 0.8. 

To summarize, we have found that in the presence of the 't Hooft interaction, treated beyound the lowest order 
SPA, several solutions to the gap equations are possible at some range of input parameters, i.e. the same values of 
G, K, A,mi lead to different sets of constituent quark masses (m„,ms) and, therefore, to different values of the quark 
condensates. A quite different scenario emerges for the hadronic vacuum, which can now be multivalued. It makes 
our result essentially different from the ones obtained in These findings must be further analysed in order to 

establish which of the extrema correspond to minima or maxima of the effective potential. This step will be done 
elsewhere in conjunction with the determination of the meson mass spectrum, as it also requires dealing with the 
terms with two powers of the meson fields in the ansatz of solutions Eq. (12^) and in the related Lagrangian (pTf) . 



IV. CONCLUDING REMARKS 



The purpose of this work has been twofold. Firstly we have developed the technique which is necessary to go 
beyound the lowest order SPA in the problem of the path integral bosonization of the 't Hooft six quark interaction. 
We have shown how the pre-exponential factor, connected with the steepest descent approach and which is responsible 
for the quantum fluctuations around the classical path, can be treated exactly, order by order, in a scheme of increasing 
number of mesonic fields, while preserving all chiral symmetry requirements. This technique is rather general and 
can be readily used in other applications. Second, we have explored with considerable detail the implications of 
taking the quantum fluctuations in account in the description of the hadronic vacuum. A very complex multivalued 
vacuum emerges at fixed values of the input parameters G, k, A and current quark masses. We encountered several 
classes of solutions. Searching in an interval of constituent quark masses from zero to ~ 3GeV, we found G, k, regions 
caracterized by one, three and four solutions. The multiple vacua may have very interesting physical consequences 
and applications. 
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The 't Hooft six quark flavor mixing interaction (iV/ = 3) is bosonized by the path integral 
method. The considered complete Lagrangian is constructed on the basis of the combined 't Hooft 
and f/(3) x L^(3) extended chiral four fermion Nambu - Jona-Lasinio interactions. The method of 
the steepest descents is used to derive the effective mesonic Lagrangian. Additionally to the Icnown 
lowest order stationary phase (SP) result of Reinhardt and Alkofer we obtain the contribution from 
the small quantum fluctuations of bosonic conflgurations around their stationary phase trajectories. 
It affects the vacuum state of hadrons at low energies: whereas without the inclusion of quantum 
fluctuations the vacuum is uniquely defined for a fixed set of the model parameters, fiuctuations give 
rise to multivalued solutions of the gap equations, marked at instances by drastic changes in the 
quark condensates. We derive the new gap equations and analyse them in comparison with known 
results. We classify the solutions according to the number of extrema they may accomodate. We 
find up to four solutions in the < m„,ms < 3 GeV region. 

12.39.Fe, 11.30.Rd, 11.30.Qc 



I. INTRODUCTION 

The global Ul{3) x Ufi{3) chiral symmetry of the QCD Lagrangian (for massless quarks) is broken by the f//i(l) 
Adler-Bell-Jackiw anomaly of the SU{3) singlet axial current q'^n'^hl- Through the study of instantons [1,2], it has 
been realized that this anomaly has physical efTects with the result that the theory contains neither a conserved U{1) 
quantum number, nor an extra Goldstone boson. Instead the effective 2Nf quark interactions arise, which are known 
as 't Hooft interactions. In the case of two flavors they are four-fermion interactions, and the resulting low-energy 
theory resembles the old Nambu - Jona-Lasinio model [3] . In the case of three flavors they are six-fermion interactions 
which are responsible for the correct description of ?y and ?y' physics, and additionally lead to the OZI-violating effects 
[4,5], 

^det = K{detqPRq + dctqPiq) (1) 

where the matrices Pr^l = (1 ± 75)/2 are projectors and determinant is over flavor indices. 

The physical degrees of freedom of QCD at low-energies are mesons. The bosonization of the effective quark 
interaction (1) by the path integral approach has been considered in [6], where the lowest order stationary phase 
approximation (SPA) has been used to estimate the leading contribution from the 't Hooft determinant. In this 
approximation the functional integral is dominated by the stationary trajectories /■g^(a;), determined by the extremum 
condition 6S{r) = of the action S{r). The lowest order SPA corresponds to the case in which the integrals 
associated with S'^S{r), for the path rg^{x) are neglected and only S'(''st) contributes to the generating functional. 
The next natural step in this scenario is to complete the semiclassical result of Reinhardt and Alkofer by including 
the contribution from the integrals associated with the second functional derivative S^S{rgi), and this is the subject 
of our paper. 

An alternative method of bosonizing the 't Hooft determinant has been reviewed in [7] . The special path integral 
representation for the quark determinant has been obtained by considering Nc as an algebraically large parameter. 
One should not forget that the 't Hooft 's determinant interactions are induced by instantons and only can be written 
in the simple determinantal form (1) in the limit of large number of colours - otherwise the many-fermion interactions 
have a more complicated structure. For our calculations it means, in particular, that the terms of the bosonized 
Lagrangian induced by the 't Hooft's determinant interaction (1) should be of order 1, corresponding to the standard 
rules of Nc counting. In this respect it is worthwhile to note that we have found that the meson vertices, induced by 
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the 't Hooft's determinant interactions in the lowest order SPA and the term from the integral associated with the 
second functional derivative PS{rgi) have the same TV^-order and should be considered on the same footing. 



II. PATH INTEGRAL BOSONIZATION 



To be definite, let us consider the theory of the quark fields in four dimensional Minkowski space, with dynamics 
defined by the Lagrangian density 

= %JL +'^det- (2) 

The first term here is the extended version of the Nambu - Jona-Lasinio (NJL) Lagrangian -CjyjL = jCq + -Cjjit' 
consisting of the free field part 

jCo = qiil^df, - m)q, (3) 
and the U{3)l x U{3)r chiral symmetric four-quark interaction 

Ant = jiiQ^aqf + {qtl^Kq?]- (4) 

We assume that quark fields have color and flavor indices running through the set i = 1,2,3; Xa are the standard 
?7(3) Gell-Mann matrices with a = 0, 1, . . . , 8. The current quark mass, rh, is a nondegenerate diagonal matrix with 
elements diag{rhu, rhd, rhs), it explicitly breaks the global chiral U{3)l x U{3)r symmetry of the jC-^jj^ Lagrangian. 
The second term in (2) is given by (1). Letting 

Sa = -qXaq, Pa = qh^Kq, S = SqAq, p = PaK (5) 

yields 

''^det = ['^^^(^ + + '^^^(^ ~ 

with determinants written in terms of the mesonic type quark bilinears. This identity is a first step to the bosonization 
of the theory with Lagrangian (2) . 

The dynamics of the system is described by the vacuum transition amplitude in the form of the path integral 

Z = j VqVqexp (^i J d^xCj . (7) 

By means of a simple trick, suggested by Reinhardt and Alkofer, it is easy to write down this amplitude as 

Z = J VqVqVaaV(j>aVrlVr^ exp [i j d'^xC'^ (8) 

with 

£' = qiiY'd;, -m-a + i-i^<p)q - ^ [{oaf + {4>a?] + r^{(Ta + GqXaq) 

+ r^((j)a - Gqi-y^Xaq) - [det(cr + + det{a - icf))] . (9) 

Eq.(8) defines the same expression as Eq.(7). To see this one has to integrate first over auxihary fields rf, rf . It 
leads to 5-functionals which can be integrated out by taking integrals over aa , and (pa , and which bring us back to the 
expression (7). From the other side, it is easy to rewrite Eq.(8) in a form appropriate to finish the bosonization, i.e., 
to calculate the integrals over quark fields and integrate out from Z the unphysical part of the auxiliary rf , scalar 
fields. Indeed, introducing new variables (T ^ (T + Gri, ^ Gr2, and after that ri —> 2ri — cr/G, r2 — > 2r2 — (/i/G 
we have 

Z = J T>(Ta'D(l)aT>qT>qexp J d'^xCq{q,q,a,(j))^ j DriQDr2Q exp j d^xCr{cJ,(j),ri,r2)^ (10) 
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where 



^, = qii-y^d^ - m - a + i'y5(j))q, (11) 
jCr = 2G [{riaf + (r2a)'] - 2(ri,(T, + r2a<^a) - I [det(ri + irj) + det(ri - ira)] . (12) 

The Fermi fields enter the action bilinearly, we can always integrate over them, because in this case we deal with the 
standard Gaussian type integral. At this stage one should also shift the scalar fields aa — > ctq + by demanding 
that the vacuum expectation values of the shifted fields vanish < 0|(Ta|0 >= 0. In other words, all tadpole graphs in 
the end should sum to zero, giving us gap equations to fix parameters Aq. Here Aq = rria — rha, with ttIq denoting 
the constituent quark masses ^. To evaluate path integrals over ri,2 one has to use the method of stationary phase, 
or, after the formal analytic continuation in the time coordinate X4 = ixo, the method of steepest descents. Let us 
consider this task in some detail. 

The Euclidean (imaginary time) version of the path integral under consideration is 



J{a,(f))= J VriaVr2aexp d^x£r{a,(t>,ri,r2) 



(13) 



This integral is hopelessly divergent even if k = 0. One should say at this point that we are not really interested in 
(13) but only in its analytic continuation. Let us suppose we analytically change £.r{cr, 4>, ri , r2) in some way such that 
we go from this situation back to the one of interest. To keep the integral convergent, we must distort the contour 
of integration into the complex plane following the standard procedure of the method of the steepest descents. This 
method gives the first term in an asymptotic expansion of J{a,(j)), vahd for h ^ 0. We lead the contour along the 
straight line which is parallel to the imaginary axis and crosses the real axis at the saddle point r^^. It is in the sense 
of this continuation that the integral J{a, </>) of (13) is to be interpreted as 

'Dria'Dr2aexp I / d'^x£.ri(^,<P,ri,r2) j ■ (14) 

Near the saddle point Tg^, 

£^ftj£^(rst) + ^^raC/3(^st)''/3 (15) 

where the saddle point, r^^, is a solution of the equations £^(ri,r2) = determining a flat spot of the surface 
C.r{rx,r2). 

r 2Gr? - (<7 + A), - f A.6c(r?rJ - r^r^) =0 , . 

\ 2Gr^ -cj)a + ^Aabcry2 = 0- ^ ^ 

This system is well-known from [6]. The totally symmetric constants, Aatc, come from the definition of the fiavor 
determinant: detr = Aabc'''""i'''r'' , and equal to 

Aabc = '^^^ijk^mnl{^a)im{^b)jn{^c)kl- (17) 

They are closely related with the L^(3) constants dabc- We use in (15) symbols f for the differences {r"^ — r^-^). To 
deal with the multitude of integrals in (14) we define a column vector f with eighteen components fa = (fijff) and 
with the matrix jC'^pir^^) being equal to 

0(r,t)=4GQ„,, Q.p=(^''^f^;f'^''' sM^aI^^^^ )■ 

V 8G"^a'":'2St "ab ^ gG^abcl -^gi / 

Eq.(14) can now be concisely written as 



^The shift by the current quark mass is needed to hit the correct vacuum state, see e.g. [8]. 
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J{a, ct>) = exp (^j d*xCr{rst)j J Vf^ exp (^G j d^xf^Q{r^^)f^ [1 + . (19) 

Our next task is to evaluate the integrals over Tq;. Before we do this, though, some comments should be made about 
what we have done so far: 

(1) The first exponential factor in Eq.(19) is not new. It has been obtained by Reinhardt and Alkofer in [6]. A bit 
of manipulation with expressions (12) and (16) leads us to the result 

'C.(rst) = I {G[{rl,0' + K.tf] " 2[('^ + A).r^",t + '^'^^2st]} • (20) 



One can try to solve Eqs.(16) looking for solutions r°g^ and r^g^ in the form of increasing powers in (Tq, 0a 

'^ist 



ha + ftiV + h^^la.o, + h^^lMc + ... (21) 



Putting these expansions in Eqs.(16) one can obtain the series of selfconsistent equations to determine the constants 
ha, /liV, and h^^^ 

2Gha -Aa- Y^abchhc = 0, (23) 
2G (^6ac - ^^acbh^ /ii'e' = ^ae, (24) 
2G (^Sac + ^^acbht^ h^c} = Sae. (25) 

The other constants can be obtained from these ones, for instance, we have 

^^abc — ~^^iah[ih[g Ag^ig, h\^^^ = - — ha^h^^^h[gAg^^g, h^^c — ~~^ha^h^^^h[gAg^^g. (26) 

As a result the effective Lagrangian (20) can be expanded in powers of meson fields. Such an expansion (up to the 
terms which are cubic in (Ja,(t>a) looks like 

Lrirst) = -^haaa - h^ab^aOb " h^aUa't>b + O(field^). (27) 

(2) Our result (19) has been based on the assumption that all eigenvalues of matrix Q are positive. It is true, for 
instance, if k = 0. It may happen, however, that some eigenvalues of Q are negative for some range of parameters 
G and n. In these cases there are no conceptual difficulties, for from the very beginning we deal with well defined 
Gaussian integrals and the integration over the corresponding fa simply does not require analytic continuation. 

We now turn to the evaluation of the path integral in Eq.(19). In order to define the measure Vfa more accurately 
let us expand fa in a Fourier series 

oo 

fa{x) = '^Cn,a^n{x), (28) 

assuming that suitable boundary conditions are imposed. The set of the real functions {^pn{x)] form an orthonormal 
and complete sequence 

/oo 
d'^X^Pn{x)^Pm{x) = 6nm, ^ V>n{x)ipn{y) = 6{x - y). (29) 



Therefore 

C 



/ Wacxp (2G / d'^xf^Q{r^{)f] = / dc„,a exp |2G V c„,„A^^Cm,^ [ = — ^ 
The normalization constant C is not important for the following. The matrix AJJ^ is equal to 



(30) 
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Kt = j d'^X¥>n{x)Qal}{x)ipm{x). (31) 

From (18) and (29) it follows that 



with 



i'aP - , Aeba \ (2) 1,(2) /„a . x I ' '^'^'^) 



'r°- ' (r°- — h \ I 



Only the matrix F^^ depends here on fields a, (j). By absorbing in C the irrelevant field independent part of 2G\^!^, 
and expanding the logarithm in the representation det ( 1 + F) = exp tr In ( 1 + F) , one can obtain finally for the integral 
in (19) 

r ( 1 °o (-1)" °° 1 

Jia, 4>) = C'e^^ Sr = d'xl £,(rst) + - ^ ^^tr [F^^ir^t)] ^ ^ra{x)^ra{x) . (34) 

•' I n=l m=l J 

The sum over m in this expression, however, is not well defined and needs to be regularized. One can regularize it by 
introducing a Gaussian cutoff M damping the contributions from the large momenta fc^ 

£ ^„(.)^„(.) = ^(0) ~ exp = ^ (35) 

This procedure does not decrease the predictability of the model, for anyway one has to regularize the quark loop 
contributions in (10). Alternatively, following ideas presented in [9], one can introduce the Ansatz: 

Sr = jd\ |/:,(rst) + ^ E ^tr [F:,(rst)] | (36) 

proposing that the undetermined dimensionless constant a will be fixed by confronting the model with experiment 
afterwards. 



III. THE GROUND STATE 



Let us study the ground state of the model under consideration, then properties of the excitations will follow 
naturally. To make further progress let us note that Eqs.(23) have non-trivial solutions for /iq, /is, /is, corresponding to 
the spontaneous breaking of chiral symmetry in the physical vacuum state with order parameters A, ^ {i = u, d, s). 
We may then use this fact to rewrite Eqs.(23) as a system of only three equations to give hi 

2Ghi - A, = '^tijkhjhk (37) 

o 

where the totally symmetric coefficients tijk are equal to zero except for tuds = 1- They are related to coefficients 
Aabc by the embedding formula SuJiaAabc^bjeck = tijk where matrices ujia, and Bq, are defined as follows 
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FIG. 1. The couplings G (left panel) and k (right panel) as functions of m„ at fixed values of nis = 572 MeV and other 
parameters in Eq.(41). These curves show the typical m„-dependence of the functions if one neglects the fluctuations term in 
Eq.(39). 

^ /^/2 ^/2 ^/2 \ ^ / 1 \ 

eai = —j= ^/l -V3 , ojia = ^\ -v^ 1 . (38) 
2V3 \^ 1 1-2 7 -2 / 

Here the index a runs a = 0, 3, 8 (for the other values of a the corresponding matrix elements are equal to zero). We 
have also ha = Saihi, and hi = ujiaha- Similar relations can be obtained for Aj and Aq. In accordance with these 
notations we will use, for instance, that /i^J-* = ujiah^j,} . 

A tadpole graphs calculation gives for the gap equations the following result 

+ ^{^- ^4«^c/4!^ = ^m,Jo(mD (39) 

where the left hand side is the contribution from (36) and the right hand side is the contribution of the quark loop 
from (10) with a regularized quadratically divergent integral Jo(to^) being defined as 

Mm'')=j^ ^e-*"'p(i,A2), p(i,A2) = l-(l + iA2)exp(-iA2). (40) 




FIG. 2. The couplings G (left panel) and k (right panel) as functions of m„ for the same (as in Fig.l) values of other flxed 
parameters and for the case when the fluctuations term in Eq.(39) is taken into account. 
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FIG. 3. Different classes of m„ solutions of the gap equations, including fluctuations, at fixed G, k values. Four branches of 
solutions, A, B, C and a — b stretch have negative k. Three branches, D, E and along c — d arm correspond to positive k. Bold 
dashes indicate that only one solution exists in this region of G, k. See further details in text. 

The second term on the left hand side of Eq.(39) is the correction resulting from the Gaussian integrals of the 
steepest descent method, comprising the effects of small fluctuations around the stationary path. If one puts for a 
moment a = in Eq.(39), and combines the result with Eqs.(37), one finds gap equations which are very similar to the 
ones obtained in [4] (see equation (2.12) therein). At fixed input parameters G,k,A of the model, the gap equations 
can be solved giving us the constituent quark masses rrii as functions of the current quark masses, rrii = mi{irij). 
Alternatively, by fixing m,, one can obtain from the gap equations the non-trivial solutions m, = mi{G, k, A). In 
particular, when = rhd, these equations can be solved for G and k, giving expressions 

~\nJ mlJl{ml)-mlJl{ml) ' ''~\nJ Jo(m2)[my2(^2) _ ^2 j2(^2)] " ^ ^ 

In Fig.l we plot the curves of G and k versus m„ keeping constant A = 0.87 GeV, = 572 MeV, = 
200 MeV, rhu = 6 MeV. One can readily see that at given values for {mu,ms) the curves yield unique values 
of (G, k), i.e. the vacuum state is well-defined in this case. 




FIG. 4. The same as in figure 3 for the rus solutions. The dashes are the constant rus corresponding to the abed curve of 
Figs. 2 and 3. 

Let us consider now the general case which we have when a > in Eq.(39). To illustrate the qualitative difference 
with the previous case we put for definiteness a = 1 and look again for the solutions G = G{mu) and k = K{mu) 
with the same set of fixed parameters. The corresponding curves are plotted in Fig. 2. If the mass m„ is sufficiently 
low that niu < ml'"*"\ (in the figures denoted by the region left the turning point c) or sufficiently high that 
^(jnax) ^ ^ ^.^ ^j^^ figures denoted by the region right to the turning point 6), then there exists again a 
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single solution with unique values of (G, k). However, there is now a region for in which rriu"' < rriu < , 
where three values of coupHngs {G, k) are possible. 

Conversely, one can study the solutions: niu = m„(G, /t), = ms{G,K at fixed values of input parameters: 
A,mj, = rhd,rhs). As starting input values for the couplings G and k in the gap equations we take the ones already 
determined along the path abed shown in Fig. 2, obtained at a constant value of the strange quark mass, rris = 572 
MeV. For a chosen value of the set {G,k) we then search for further solutions {mu,ms) of the Eq.(39), displaying 
results in Fig. 3 for niu and in Fig. 4 for m^, correspondingly. The dashed curves are the repetition of the solutions 
encountered in Fig. 2. The bold dashes in Fig. 3 indicate that we only find this one solution at fixed values of G, k. 
Combining the information of Figs. 3 and 4, one sees that one has up to four solutions at fixed G, k. Indeed, travelling 
along the original path abed one observes the following: the branch ab is accompanied by three other branches, marked 
as A, B and C which belong to the same class of solutions. One sees that these solutions have negative k values. From 
the turning point b until the maximum value of G (and corresponding k) (bold dashes) we have no other solutions to 
the gap equations, as already stated. From this maximum G value to the turning point c and further along the ed 
arm we encounter further two branches, denoted by D and E to the solutions of Eqs.(39) with same G, k. They are 
positive K solutions. 

This very rich structure of the vacuum solutions implies the possibility of having several different values of the quark 
condensates for the same G, k parameters, embracing also the possibility which has been considered in connection 
with generahzed chiral perturbation theory (see e.g. section 4 in [10]). We give here only a few examples. At 
G = 4.54 GeV"^,K = 153.04 GeV"^ we have three solutions. The solution m„ = 346 MeV, = 572 MeV, on 
the ed arm has the quark condensates < uu >i/^= —236.8 MeV, < ss >i/^= —183.5 MeV and the ratio R = 

1 /3 

{< ss > / < uu >) ' = 0.775; the second solution, on the E branch with m„ = 535 MeV,ms = 732 MeV has 
condensates < uu >'^/^= -249 MeV, < ss >'^/^= -183 MeV and the ratio R = 0.738; the third solution, located at 
the D branch with niu = 46 McV,ms = 418 MeV has condensates < uu >i/3= _i3i MeV, < ss >V3= _i72 MeV 
and R = 1.313. We chose the next example at G = 8.126 GeV^^, k = —544.81 GeV^^, where there are four solutions. 
The solution m„ = 316 MeV,ms = 572 MeV, on the ab arm has the quark condensates < uu >^/^= —233 MeV, 

< ss >i/3= _184 MeV and R = 0.787; the second solution, on the B branch with = 624 MeV,ms = 205 MeV 
has condensates < uu >i/3= -249 MeV, < ss >i/2= -56.7 MeV and R = 0.227; the third solution, located at the 
A branch with m„ = 421 MeV, iris = 1-84 GeV has condensates < uu >i/3= -14.4 MeV, < ss >i/3= -86 MeV 
and R = 0.353. The fourth solution, on the C branch, with to„ = 1.394 GeV,TOs = 1.594 GeV has condensates 

< uu >^/3= -230 MeV, < ss >^/^= -120 MeV and the ratio R = 0.524. As a final example we take the solutions 
at G = 11.96 GeV"^ ~ Gmax,K = 371.491 GeV"^ where the branches D and E emerge and are very close to 
each other with m„ = 1.56 GeV,ms = 1.72 GeV. The corresponding condensates are < uu >i/^= —224 MeV, 

< ss >i/3= -105 MeV and R = 0.467. The other solution is at the path be with m« = 290 MeV, = 572 MeV 
with condensates < uu >V3= -229 MeV, < ss >V3= _184 MeV and i? = 0.8. 

To summarize, we have found that in the presence of the 't Hooft interaction, treated beyound the lowest order 
SPA, several solutions to the gap equations are possible at some range of input parameters, i.e. the same values of 
G,K,A,rhi lead to different sets of constituent quark masses [niu, iris) and, therefore, to different values of the quark 
condensates. A quite different scenario emerges for the hadronic vacuum, which can now be multivalued. It makes 
our result essentially different from the ones obtained in [6,4]. These findings must be further analysed in order to 
establish which of the extrema correspond to minima or maxima of the effective potential. This step will be done 
elsewhere in conjunction with the determination of the meson mass spectrum, as it also requires deahng with the 
terms with two powers of the meson fields in the ansatz of solutions Eq.(21) and in the related Lagrangian (27). 

IV. CONCLUDING REMARKS 

The purpose of this work has been twofold. Firstly we have developed the technique which is necessary to go 
beyound the lowest order SPA in the problem of the path integral bosonization of the 't Hooft six quark interaction. 
We have shown how the pre-exponential factor, connected with the steepest descent approach and which is responsible 
for the quantum fluctuations around the classical path, can be treated exactly, order by order, in a scheme of increasing 
number of mesonic fields, while preserving all chiral symmetry requirements. This technique is rather general and 
can be readily used in other applications. Second, we have explored with considerable detail the implications of 
taking the quantum fluctuations in account in the description of the hadronic vacuum. A very complex multivalued 
vacuum emerges at fixed values of the input parameters G, k, A and current quark masses. We encountered several 
classes of solutions. Searching in an interval of constituent quark masses from zero to ~ 3GeV, we found G, /t, regions 
caracterized by one, three and four solutions. The multiple vacua may have very interesting physical consequences 
and applications. 
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